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Abstract 

Two-dimensional fluid simulations of interchange turbulence for geometry and parameters relevant for 
the scrape-off layer of confined plasmas are presented. We observe bursty ejection of particles and heat 
from the bulk plasma in the form of blobs. These structures propagate far into the scrape-off layer where 
they are lost due to transport along open magnetic field lines. From single-point recordings it is shown 
that the blobs have asymmetric conditional wave forms and lead to positively skewed and flat probability 
distribution functions. The radial propagation velocity may reach one tenth of the sound speed. These 
results are in excellent agreement with recent experimental measurements. 

PACS numbers: 52.25.Gj, 52.35.Ra, 52.65.Kj 
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Recently, several experimental investigations have revealed a strongly intermittent nature of 
particle and heat transport in the scrape-off layer (SOL) of magnetized plasmas Jlll^lsl. There 
are strong indications that this is caused by localized structures in the form of plasma "blobs" 
propagating radially far into the SOL. It has been suggested that this is due to the dipolar vorticity 
field caused by vertical guiding center motions in plasmas in non-uniform magnetic fields fl], 
reflecting the compressibility of the diamagnetic current. An outstanding challenge is to give a 
self-consistent description of the emergence and evolution of such structures, which also capture 
their statistical properties. Here an attempt towards this goal is presented, yielding favorable 
agreement with experimental measurements. This is achieved by focusing on the collective two- 
dimensional dynamics perpendicular to the magnetic field while using a simplified description of 
particle and heat losses to limiters or end plates along open field lines. 

Some of the most prominent features of experimental single-point measurements are the pres- 
ence of asymmetric conditional wave forms as well as strongly skewed and flat probability distri- 
bution functions (PDF's) of the density and temperature signals fly,!^]. Simple interpretations as 
well as advanced imaging techniques give a picture of field-aligned blobs or filaments propagating 
out of the bulk plasma with radial velocities up to one tenth of the sound speed. These highly 
non-linear thermal structures have amplitudes which significantly exceed the back- ground levels. 
The associated intermittent transport may have severe consequences for magnetic confinement ex- 
periments by producing large heat bursts on plasma facing components. A change of sign in the 
asymmetry of the fluctuation time series PDF close to the last closed flux surface (LCFS) indicates 
that the structures are generated by a "flapping" of the edge pressure gradient, ejecting blobs of 
excess particles and heat out of the bulk plasma iQ]. The blob structures thus seems to be generated 
close to the LCFS and subsequently propagate far into the SOL, where they are subject to strong 
damping due to parallel losses. This separation of driving and damping regions in configuration 
space has been discarded in several previous studies of SOL turbulence o,KJ]. 

In this Letter we present a novel model for interchange turbulence in slab geometry and nu- 
merical solutions in qualitative agreement with experimental measurements. The model geometry 
comprises distinct production and loss regions, corresponding to the edge and SOL of magnetized 
plasmas. The separation of these two regions defines an effective LCFS, though we do not include 
magnetic shear in our model. In the edge region, strong pressure gradients maintain a state of 
turbulent convection. It is demonstrated that a self-regulation mechanism involving differential 
rotation leads to a repetitive expulsion of hot plasma into the SOL, resulting in asymmetric con- 
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ditional wave forms, strongly non-Gaussian probability distributions and significant cross-field 
transport by localized objects. 

Assuming cold ions and neglecting electron inertia effects, a three-field model may be derived 
for quasi-neutral electrostatic perturbations of the full particle density n, electric potential § and 
electron temperature T. Using the Bohm normalization and slab coordinates with z along the 
magnetic field we obtain 021], 

3 +?xVfvk-c(p) = vnV 2 Sl-a Q £l, 



dt 
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where time is normalized by the ion gyration period, 1 / co c ,, and spatial scales by the hybrid gy- 
ration radius, p s = c s /co c ,-. The particle density n and temperature T are normalized to fixed 
characteristic values at the outer wall. We further define the two-dimensional advective derivative, 
the magnetic field curvature operator and the toroidal magnetic field, respectively, by 



dt dt B dy \+£ + tpc 

The vorticity is given by £1 = V^(|), the inverse aspect ratio £ = ci/Rq and C, = p s /Ro where a and 
Rq are the minor and major radius of the device, respectively. The terms on the right hand side 
of the model equations describe external sources 5, parallel losses along open field lines through 
the damping rates a, and collisional diffusion with coefficients v. The geometry and boundary 
conditions are sketched in Fig.[T] 

In the absence of external forcing and dissipative processes the model equations non-linearly 
conserves the global energy 

E = Jdx I(v ± ^) 2 + ^r 

where the integral extends over the whole plasma layer under consideration. Thus, the curvature 
terms due to magnetic field inhomogeneity correctly yield a conservative energy transfer from the 
confined heat to the convective motions. We further define the kinetic energy of the fluctuating 
and mean components of the flows, 

K = Jdx 1 -(v A _^), U = Jdx~vl (1) 
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where the zero index denotes an average over the periodic direction y and the spatial fluctuation 
about this mean is indicated by a tilde. The linearly damped mean flows, vo = dtyo/dx, does not 
yield any radial convective transport and hence form a benign path for fluctuation energy. The 
energy transfer rates from thermal energy to the fluctuating motions, and from the fluctuating to 
the mean flows, are given respectively by 

F p = j dxnTc(<$>), F v = J dxv x v y ^. (2) 



momentum up the 

islflH. 



Note that F p is essentially a measure of the domain integrated convective thermal energy trans- 
port, while F v shows that structures tilted such as to transport positive poloida 
gradient of a sheared flow will sustain the flow against collisional dissipation 

In the following we present results from numerical simulations of the interchange model using 
parameters relevant for SOL plasmas. The dimensions of the simulation domain is L x = 2L y = 400 
and the LCFS is located at x LCFS = 100. The parameters used for the simulation presented here are 
£ = 0.25, C, = 10~ 3 , and the collisional diffusion v = 5 x 10 -3 is taken to be the same for all 
fields. The parallel loss rate of temperature is assumed to be five times stronger than that on the 
density and vorticity, o„ = Oq, = Gr/5 = 10~ 3 /2ftg, since primarily hot electrons are lost through 
the end sheaths [11]. The damping rates o„ and Oq correspond to losses along one connection 
length InRoq with the acoustic speed c s , where q = 3 is the safety factor at the edge. Finally, the 
radial line-integral of the sources S n and St equals 0.2, and the shape of the sources and parallel 
loss coefficients shown in Fig. [T] are given by 8 = 16 and £, = 2. For the numerical solution we 
have employed an Arakawa scheme for the advective non-linearities and a third order stiffly stable 
scheme for the time integration [Q]. The spatial resolution is 512 and 256 grid points in the radial 
and poloidal directions, respectively, and the total time span of the simulation is 4 x 10 6 . 

In Fig. [2] we show the typical evolution of the particle confinement P and heat confinement H 
in the edge and SOL regions, defined by 

f-^LCFS fL x 



/"*LCFS rL x 

/ dxn (x,t), P S ol= / dxn (x,t 

J0 JXir^s 



*LCFS 

and similarly for the heat confinement H. From the figure we observe that plasma and heat gradu- 
ally builds up in the edge at the same time as it is decaying in the SOL region. This is repetitively 
interrupted by rapid changes in which plasma and heat is lost from the edge to the SOL. More than 
20% of the edge plasma may be lost during individual bursts. Also note from the figure that the 
normalized heat confinement is much less than the particle confinement due to the larger loss rate 



in the SOL region of the former. Further insight is revealed by Fig. |3] which shows the evolution 
of the kinetic energy contained by the mean and fluctuating motions, confer Eq. (U]), as well as the 
collective energy transfer terms defined in Eq. ©. From the figure we observe that the convective 
energy and thermal transport appears as bursts duringwhich particles and heat are lost from the 
edge into the SOL region. As discussed in Refs. [^,0, [ill this global dynamics is caused by a 
self -regulation mechanism in which kinetic energy is transfered from the fluctuating to the mean 
components of the flows, and subsequently damped by collisional dissipation. As the thermal con- 
finement is allowed to vary, this results in characteristic sawtooth oscillations |8]. Note, however, 
the clear demonstration in Figs. El and |3] that thermal energy is ejected in a bursty manner from the 
edge and into the SOL region, where it is eventually lost by transport along open field lines. 

The time-averaged profiles of the plasma density and the mean poloidal flows are shown in 
Fig. IU Parallel losses in the SOL region result in average profiles peaked inside the LCFS, and 
weakly decaying throughout the SOL. Also shown in Fig. |4] are typical instantaneous profiles 
during a quiet period (/ qu iet) and during a burst (/burst)- We observe significant deviations from the 
average profiles, with a more peaked density profile in the edge region during quiet phases. During 
bursts there is a substantial increase of the density profile in the SOL region due to the convective 
plasma transport. The temperature profiles have a similar structure but with lower amplitudes due 
to the larger loss rate. The self-sustained poloidal flow profiles are strongly sheared in the edge 
region, and have larger amplitudes during the the strong fluctuation period. This figure clearly 
indicates that fluctuations are driven in the strong pressure gradients in the edge region. 

The statistics of single-point recordings at different radial positions Pi indicated in Fig. incom- 
pletely agree with experimental measurements. In Fig. 13 we present the probability distribution 
functions of the density signals taken from a long-run simulation containing more than a hundred 
burst events. The flat and strongly skewed distributions indicate the high probability of large pos- 
itive fluctuations corresponding to blobs of excess plasma. The skewness and flatness factors take 
values up to 3 and 15, respectively, except for the outermost radial position (Pj) where very few 
structures arrive. At the other points the PDF's have similar structure with a pronounced expo- 
nential tail towards large values, a characteristic feature of turbulent convection in the presence of 
sheared flows flU 

The conditionally averaged temporal wave forms calculated from the same signals, using the 
trigger condition n > 4n ims at each individual point, are presented in Fig. [6] An asymmetric wave 
form with a sharp rise and a relatively slow decay, as observed in experimental measurements, is 
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clearly seen yj,|2j,|3|]. The maximum density excursions significantly exceed the background level, 
and decay rapidly as the structures propagate through the SOL. The number of realizations for the 
conditional averaging decreases gradually from 125 at the innermost probe to 2 at the outermost 
one. Using a negative amplitude for the conditional averaging results in very few realizations, 
again showing the presence of blob-like structures. From two-dimensional animations and simple 
statistical correlation measurements we further find that the radial propagation velocity of these 
structures is typically about one tenth of the sound speed but with a large statistical variance. This 
is in excellent agreement with experimental measurements JllQ,Q]. 

In Fig. 13 we show the spatial structure of the density, temperature, vorticity and electrostatic 
potential during a quiet period and during a burst. These correspond to the same times as the 
instantaneous profiles shown in Fig. HJ hi the quiet period there are only weak spatial fluctua- 
tions with the plasma and heat well confined within the LCFS. During bursts we observe strong 
structures in all fields, which have propagated far into the SOL. Notice again the much weaker 
perturbation in the temperature field as compared to the density due to the difference in the par- 
allel loss rates. Moreover, while blob-like structures are observed for the density and temperature 
fields, the vorticity displays a roughly dipolar structure as expected from theory and experimen- 
tally measured yj,|4|]. 

In this Letter we have proposed a new model for interchange turbulence and demonstrated that 
its numerical solutions are in good agreement with that reported from experimental investigations 
of SOL turbulence and transient transport events An important feature of our model 

is the spatial separation between forcing and damping regions which has not been accounted for 
in previous studies |Q, E]. While our model does not in detail describe the transition region be- 
tween closed and open magnetic field lines, its solution clearly demonstrates a sound mechanism 
for the origin and nature of intermittent transport in the SOL of magnetized plasmas. Our work 

verifies the present experimental working hypothesis in terms of field-aligned blob-like structures 

fj 

propagating far into the scrape-off layer, which was recently questioned by Sarazin et al. |5]. The 
blob-like transport results in strongly skewed and flat PDF's and asymmetric conditional wave 
forms. This is caused by a "flapping" of the edge pressure profile, ejecting blobs of excess par- 
ticles and heat out of the bulk plasma. The excellent qualitative agreement between simulation 
data and experimental observations gives strong confidence that the two-dimensional structure in 
potential, vorticity, density and temperature of the blobs reveals the actual spatial shape of these 
blobs as occurring in experiments. 
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FIG. 1: Geometry of the simulation domain showing the forcing region to the left, corresponding to the 
edge plasma, and the parallel loss region to the right, corresponding to the scrape-off layer. Data time series 
are collected at the probe positions p. 

FIG. 2: Evolution of particle confinement P and the heat confinement H in the edge and scrape-off layer 
regions, showing sawtooth oscillations. 

FIG. 3: Evolution of the kinetic energies and the collective energy transfer terms, showing bursty behavior 
in the fluctuation integrals. 

FIG. 4: Time-averaged profiles of plasma particle density n and mean poloidal flow vo, and typical profiles 
during a fluctuation burst and during a quite period. 

FIG. 5: Probability distribution functions of density measured at seven different radial positions P, as shown 
in Fig. [2 The vertical axis shows count numbers. 

FIG. 6: Conditionally averaged wave forms of the density measured at seven different radial positions P, as 
shown in Fig.[0 using the condition n(xp i ) > 4w rms . 

FIG. 7: Typical spatial structure of density, temperature, vorticity and electric potential during a quite period 
to the left and during a burst to the right. 
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